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1. Introduction 

Quantum groups first appeared in the study of quantum integrable systems. Since 
then, they have proven useful in many applications, including among others 
conformal field theory, statistical mechanics, representation theory and the theory of 
hypergeometric functions, along with exhibiting a rich internal structure. Quantum 
groups have appeared seminally in the mathematical physics literature in connection 
with topological quantum field theory [I]. More recently, the quantum group su g (2) 
was used to construct "loop" [3] and "spin foam" [3J models of quantum gravity with 
a positive cosmological constant. 

The 6j-symbol (or Racah-Wigncr coefficients and henceforth the 6j or the q-6j 
for the classical and q-dcformcd versions respectively) has been studied extensively. 
It first appeared in work on g-hypergeometric functions [4]. Later, it was found 
to play a central role in the representation theory of su 9 (2) [5j. It is known to 
satisfy some recurrence relations [HE], including a particular linear, three term, single 
argument one [7J [5] ■ This recurrence has been used to analyze the asymptotics of 
the classical [9] [10] and quantum [11] 6j. 

However, this q-6j and the need to know its numerical values for a large number 
of arguments also appear in other segments of the literature [31 [T2] , less aware of these 
properties. An explicit formula for the q-6j is well known and involves a number of 
arithmetic operations that is linear in its arguments (see Appendix). In applications 
where a large number of 6j-symbols is needed at once, e.g., for all values of one 
argument with others fixed as is the case in [3] , the total number of operations becomes 
quadratic in the arguments. The above mentioned recurrence relation can be used to 
greatly increase the efficiency of the calculation by reducing the total operation count 
to be linear in the arguments. The main goal of this paper is to show that, moreover, 
this recurrence can be realized as a tridiagonal, symmetric eigenvalue problem, a 
property not shared by most recurrence relations, especially since the equivalence is 
established using only rational operations. 

A significant advantage of the eigenvalue formulation is the ability to make use 
of readily available, robust linear algebra packages, such as LAPACK [13], which 
automatically take care of the important issues of numerical accuracy and stability. 
When q = 1 or when q is a primitive root of unity, the relevant inner product becomes 
cither positive- or negative-definite and standard, specialized numerical methods can 
be exploited to increase the efficiency of the calculation even further. A secondary 
goal of this paper is to concisely collect all the relevant information needed to readily 
implement such an efficient q-6j numerical evaluation algorithm without intimate 
familiarity with the literature on quantum groups or q-hypergeometric functions. 

Sections [2] and [3j which can be skipped by those familiar with the mathematical 
literature on su g (2) recoupling theory, introduce the basic notions of the spin 
network formalism [TH [TU [16], define the Kauffman-Lins convention for the 6j- 
symbol and summarize basic diagrammatic identities needed for Section [4[ where 
the recurrence relation is realized as an eigenvalue problem. This is accomplished as a 
byproduct of an alternative proof of the recurrence itself that generalizes the classical 
argument from the Appendix of |10j . The diagrammatic spin network formalism 
makes all intermediate calculations easy to check and reproduce. The Appendix 
conveniently summarizes all formulas needed for a direct computer implementation 
of the recurrence-based evaluation of the q-6j, including the connection between the 
Kauffman-Lins and Racah-Wigner notational convention, which is traditionally used 
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in physics. 

2. Spin networks 

In a variety of physical and mathematical applications, one often encounters tensor 
contraction expressions of the form 



where T, A, B, . . . , Z are invariant tensors, with each index transforming under 
a representation of a group or an algebra. The application at hand usually calls 
for evaluating T, or at least simplifying it. An extensive literature on this subject 
exists for the classical group SU(2) or its Lie algebra, a subject known as angular 
momentum recoupling |171 118] . It is well known that such tensor contractions can be 
very efficiently expressed, manipulated, and simplified using diagrams known as spin 
networks |14j . Extensions of these techniques [151 116] are also known for the quantum 
(or q-deformed, since they depend on an arbitrary complex number q ^ 0) analogs, 
the quantum group su 9 (2) or U q (su(2)). The basics of this diagrammatic formalism, 
as needed for the derivation of the recurrence relation, are given in this and next 
sections. All relevant formulas, including explicit spin network evaluations in terms 
of quantum integers are listed in the Appendix. 

Single spin networks are edge-labeled gra Phi, where each vertex has valence 
cither 1 or 3. General spin networks are formal linear combinations of single spin 
networks. Edges attached to univalent vertices are called free. Spin networks without 
free edges are called closed. Conventionally, the labels are either integers (spins) of 
half-integers (twice- spins), which correspond to irreducible representations of su g (2). 
Reference [15] labels all spin networks with twice-spins. Unless otherwise indicated, 
all conventions in this paper follow |15j . Two spin networks may be equal even if not 
represented with identical labeled graphs. A complete description of these identities 
are given in [15] and j!6j : their study constitutes spin network recoupling and is what 
allows us to equat^§| spin networks with su 9 (2)-invariant tensors and their contractions. 

In this correspondence, each index of a tensor, transforming under an irreducible 
representation, corresponds to a spin network edge, labeled by the same representation 
(free indices correspond to free edges). In particular, a closed spin network corresponds 
to a complex number. Spin networks form a graded algebra over C (as do tensors). 
The grading is given by the number of free edges (free indices) and the product is 
diagrammatic juxtaposition (tensor product). 

3. Diagrammatic identities 

The spin networks with n free edges with fixed labels (n-valent spin networks) form 
a linear space with a natural bilinear form (or inner product). Suppose that the 
free edges are ordered in some canonical way, then, given two spin networks, we can 
reflect one of them in a mirror and connect the free edges in order. The value of 
the resulting closed spin network defines the bilinear form, which is symmetric and 

X Spin networks are actually ribbon graphs, but since all diagrams in this paper are planar, the 
ribbon structure can be added through blackboard framing. 

§ A down-to-earth guide to this correspondence, for the classical q = 1 case, can be found in 
Appendix A of I19| . Complete details with proofs can be found in I16| . 
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non-degenerate [T5]. We use the bra-ket notation for this inner product (s'|s), where 
s and s 1 are two spin networks. We also let \s) stand for s and (s'\ for the reflection 
of s' . The existence of an inner product allows the following identities, whose proofs 
can be found in [TS]- For each identity, the corresponding well known fact of SU(2) 
representation theory is given. 

The space of 2-valent spin networks, with ends labeled a and &, is 1-dimensional 
if a = b, and O-dimensional otherwise. For non-trivial dimension, the single edge gives 
a complete basis and therefore the bubble identity: 

This identity the diagrammatic analog of Schur's lemma for intcrtwiners between 
irreducible representations. 

The space of 3-valent spin networks, with ends labeled a, b and c, is also 1- 
dimensional if the triangle inequalities (|A.6|) and parity constraints (|A.7[) are satisfied, 
and O-dimensional otherwise, if q is generic. When q is a primitive root of unity, 
the dimension also vanishes whenever the further r-boundedness constraint (| A.8|) is 
violated. In the case of nontrivial dimension, the canonical trivalent vertex gives a 
complete basis and therefore the vertex collapse identity: 



a 




The normalization of the vertex, the value of the ^-network, is evaluated in 
Equation (|A.5|) . This identity is the diagrammatic analog of the uniqueness (up to 
normalization) of the Clebsch-Gordan intertwiner. 

Now, consider the space of 4-valent networks with free edges labeled a, b, c and 
d. There are two natural bases, the vertical (Z| and the horizontal \j): 




The admissible ranges for j and I, the dimension n of this space, and the conditions 
on (a, b, c, d) under which n > are given by Equations ([ATT6| through (|AT27|) . The 
transition matrix between the two bases is given by the so-called Tet-network: 



Tet(a, 6, c, d; j, I) = 
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The coefficients expressing the vertical basis in terms of the horizontal one define the 
6j-symbol, which can be expressed in terms of the Tet-network: 




Note the subscript KL for Kauffman-Lins, since this 6j-symbol is defined with respect 
to the conventions of [15] . The relation to the classical Racah-Wigner 6j-symbol used 
in the physics literature is given explicitly in Equation IA. 281 



4. Recurrence relation for the Tet-network 



The identities given in the previous section allow an alternative, elementary derivation 
of the three-term recurrence relation for the Tet-network, distinct from the standard 
one. The standard derivation is given in [7] and another is possible using the general 
theory of recurrences for g-hypergeometric functions [B], but neither directly yields 
the symmetric eigenvalue problem form. 

It is easy to check, using the bubble identity, that both the vertical and horizontal 
bases are orthogonal and that they are normalized as 



b a 




Curiously, when these normalizations are fully expanded using formulas from the 
Appendix, they take the form {—) a P/Q, where P and Q are products of positive 
quantum integers. In both cases, er = (a + b + c + d) /2, is an integer independent of j 
or I. When q = 1 or when q is a primitive root of unity, positive quantum integers arc 
positive real numbers. Hence the above inner product is real and either positive- or 
negative-definite. On the other hand, for arbitrary complex q, the normalizations (|9]) 
and (|10j) can be essentially arbitrary complex numbers. 

If we can find a linear operator L that is diagonal in one basis, but not in the 
other, then we can obtain as matrix elements of the diagonalizing transformation. 
Furthermore, if the non-diagonal form of L is tridiagonal, then the linear equations 
defining (j\l) reduce to a three-term recurrence relation. 
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Wc can construct such an operator by generalizing the argument for the classical 
case, found in the Appendix of [TO]. For brevity of notation, wc introduce a special 
modified version of the trivalent vertex: 



a I a a 



(11) 



The unlabeled edge implicitly carries twice-spin 2 and the bold dot indicates the 
multiplicative factor of [a]. Using it, we can define a symmetric operator L. Its 
diagrammatic representation and its matrix elements are given below. 

The operator L is diagonal in the |Z)-basis and its matrix elements Lw = (l\L\V) 

are 



L,, 



v — 



a m a 

I c 

\(a,b,l)(l\l'), 




with 



[ q-b+H [ -q+b+H _ [ q+b-H r 

X(a,b,l) = L 2 J L 2 J " 



a+b+l 



2] 



[2] 



where we have evaluated Tet(a, a, 6, b; I, 2) as 



(12) 
(13) 

(14) 




A(et, b, I). 



(15) 



This result may be obtained directly from Equation (|A.9[) . where the sum reduces to 
two terms, or from more fundamental considerations [20] . In the limit, q — > 1, the 
eigenvalues simplify to A(a, b, I) = + 2) — a(a + 2) — 6(6 + 2)], which shows 

that the operator L is closely related to the "square of angular momentum" in 
quantum mechanics, which was used to obtain the classical version of this recurrence 
relation [TO] . 

On the other hand, in the \j) basis, the operator L is not diagonal and the matrix 
elements Ljf = (J\L\j'}, making use of the vertex collapse identity, are 




(16) 



with the special case Loo = 0. Fortunately, though Ljji is not diagonal, it is 
tridiagonal. This property is a consequence of the conditions enforced at the central 
vertex in both Tet-networks above: the triangle inequality, \j < 2, and the parity 
constraint, which forces admissible values of j to change by 2. If these conditions are 
violated, the matrix element Ljji vanishes. 
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The diagonal elements Ljj can be evaluated using (fT5")) . For the off-diagonal 
elements Lj +2 .j = Lj.j+2 we also need 




1 



a + d — j 



j + 2 




(17) 



which can be obtained in the same way as (| X5|) . Finally, we need the identities 

2 2 

[J + 2] 

L [2M 




(18) 



The |j)-basis matrix elements can now be expressed as (again, recall the special case 
Loo = 0) 



Ljj = -O'li) 
lj,j+2 = (i+2|i+2) 



.. [2]A(a,j ) d)A(6,j,c) 



[?][?+ 2] 
a + d — j 



b + c - j 



(19) 
(20) 



The transition matrix elements can now be obtained by solving an eigenvalue 
problem in the |j)-basis: 



(i'li') 



(j'Ij') 



(21) 
(22) 
(23) 



where A; = \(a, b, I). Since Ljj> is tridiagonal, we obtain a three-term recurrence 
relation for the (j\l) transition coefficients. Expanding the expression for Ljj>, we find 
the following general form of the recurrence relation: 



(j - 2|i - 2) 



(i - 2|Z> + 



33 



\Q\~3) 



J 3J+2 



(j + 2\j + 2) 



(j + 2\l) = 0, (24) 



with the provision that Ljy vanishes whenever either of the indices fall outside the 
admissible range or j = j' = 0. Finally, the transition coefficients are uniquely 
determined (up to sign) by requiring the normalization condition 



T <mm =m 



(25) 



Practically, it is more convenient to recover the correct normalization for all j and 
fixed I, or vice versa, by requiring (j\l) to agree with (|A.9[) for j = j, cf. (|A.16[) . where 
the sum reduces to a single term. 
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Once the Tet-network has been evaluated recursively, the 6j-symbol can be 
obtained from Equation ([7]). Alternatively, a linear, three-term recurrence relation 
directly for the 6j-symbol follows from (|24[) and the linear, two-term recurrence 
relations for the bubble and ^-networks, obvious from (|A.4I) and (|A.5[) . However, 
because of the additional normalization factors in Equation ([7]) , this direct recurrence 
relation cannot be cast in the form of a symmetric eigenvalue problem like (|23|) using 
rational operations alone. 
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Appendix A. Formulas 

For a complex number q ^ and an integer n the corresponding quantum integer is 
defined as 

r q n - q- n 



q-q 1 



(A.1) 



In the limit q —y 1, we recover the regular integers, [n] — > n. When q — exp(i7r/r), for 
some integer r > 1, it is a primitive root of unity and the definition reduces to 

= sin^r/r) 
sin(7r/r) 

This expression is clearly real and positive in the range < n < r. Quantum factorials 
are direct analogs of classical factorials: 

[0]! = 1, [n]! = [l][2] .-[»]. (A.3) 

Next, we give the evaluations of some spin networks needed in the paper. They 
are reproduced from Ch. 9 of [15j . The bubble diagram evaluates to 



j= (-y\j + i] (A.4) 

whenever it is non- vanishing. For generic q, it vanishes if j < and if q is a primitive 
root of unity then it also vanishes when j > r — 2. The 8-network evaluates to 

(-)*[* + l]\[ S -a]\[ S -b]\[s-c]\ 

0{a,b,c) = , (A.5 

[a]![6]![c]! 

with s = (a + b + c)/2, whenever the twice-spins (a, b, c) are admissible and vanishes 
otherwise. Admissibility consists of the following criteria (besides the obvious a, &, c > 
0): 

{a<b+c 
b < c + a , (A.6) 
c<a + b 

parity a + b + c = Q (mod 2). (A. 7) 
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When q is a primitive root of unity, further constraints needs to be satisfied: 

f a,b,c<r-2 
r-boundedness < . (A. 8) 

|a + 6 + c<2r-4 

The tetrahedral- or Tet- network evaluates to 

^^**°-l^ aP-K-qi - (A ' 9) 

where the summation is over the range m < S < M and 

H = JJ[6 a -a t ]!, £! = [a]![6]![c]![rf]![j]![/]!, (A.fO) 

oi = (o + d + i)/2, 61 = (6 + d + j + i)/2, (A.fl) 

a 2 = (6 + c + j)/2, 6 2 = (a + c + j + l)/2, (A.12) 

as = (o + 6 + /)/2, 6 3 = (a + & + c + rf)/2, (A.13) 

a 4 = (c + d + Z)/2, m = maxja,}, (A. 14) 

M = min{6j}. (A.15) 

The indices 1 and j fully span the defined ranges. Each of the triples of twice-spins 
(a, 6,/), (c, d, Z), (a, d, j) and (c, 6, j) must be admissible, otherwise the Tet-network 
vanishes. Then, due to parity constraints, a t , 6,, m, M, and S are always integers. If 
the twice-spins (a, 6, c, d) are fixed, the admissibility conditions for generic q enforce 
the ranges of j < j < j and I < I < I to 

j = max{|a — d|, |6 — c|}, j = min{a + d, 6 + c}, (A. 16) 

Z = max{|a — 6|, |c — d|}, I = min{o + 6, c + d}, (A. 17) 

with 

j = a + 6 = c + d (mod 2), (A.18) 

Z = a + d = fe + c (mod 2). (A.19) 

The number of admissible values is the same for j and I and is equal to n = max{0, n}, 
where 

h = min{m, s — M} + 1, m = min{a, 6, c, d}, (A. 20) 

s = (a + 6 + c + d)/2, M = max{a, b, c, d}. (A. 21) 

This number n is also the dimension of the space of 4-valent spin networks with fixed 
twice-spins (a, b, c, d) labeling the free edges. This dimension is non- vanishing, n > 0, 
precisely when the twice-spins satisfy the conditions 

a + b + c + d < 2max{a, b,c, d}, (A. 22) 

a + fe + c + d = (mod 2). (A.23) 

When q is a primitive root of unity, the admissible ranges shrink to I < j < j r 
and l< I <ln where 

j r = min {j, r - 2, 2r - 4 - max{a + d, 6 + c}} , (A. 24) 

l r = min {7, r - 2, 2r - 4 - max{a + 6, c + d}} . (A. 25) 
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The number of admissible values in each range is thus restricted to n — max{0,n r }, 
where 

h r = min {fi, r — 1 — max{M, s — m}} . (A. 26) 

The condition n > requires (|A.22|) . (|A.23[) and 

a + b + c + d<2 min{a, 6, c, d} + 2r - 4. (A. 27) 

The above admissibility criteria are well known. However, the consequent explicit 
expressions for the constraints on (a, b, c, d), the bounds on j and I, and the dimension 
n are not easily found in the literature. 

In the classical q = 1 case, the Kauffman-Lins version of the 6j-symbol ([7]) differs 
from the Racah-Wigner convention used in the physics literature, which preserves the 
symmetries of the underlying Tet-network. The two 6j-symbols are related through 
the formula 

fji/2 j 2 /2 j 3 /2l Tct(Ji, J 2 ,ji,h; J3J3) 

\Ji/2 J 2 /2 Js/2j RW y/\e(J 1 ,J2,h)0(ji,32,33)6(Ji,32, WJ2J1, J 3 )| ' 

(A.28) 

which can be obtained by comparing the explicit expressions (|A.9j) and (6.3.7) of [17] . 
Note that the argument of the absolute value under the square root has sign (— y 3 ~ J s . 
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